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Int roductlon:  We  shall  assume  throughout  the  following  model 

0:Y  = XfJ+e 

where  Y is  an  n x 1 vector  of  observations,  X is  an 
n X p known  matrix  of  rank  r,  0 is  a p x 1 vector  of 
unknown  parameters,  and  e is  distributed  N(0,  a I). 

Definition  1.  A linear  parametric  function  Y = c^p,  where 
c is  a vector  of  known  constants,  is  said  to  be  an  estimable 
linear  parametric  function  if  there  exists  an  unbiased  linear 

/ 9 

estimator  a Y,  i.e.,  such  that  E(a  Y)  = Y. 

The  estlmability  of  'V  solely  depends  on  the  design 
matrix  X as  the  following  known  tests  for  estlmability  of 
Y show; 

(i)  Y is  estimable  if  and  only  if  c is  in  the  row 

space  of  X,  i.e.,  if  and  only  if  there  exists  a 

/ / 

vector  t such  that  c = t X. 

(ii)  Y is  estimable  if  and  only  if  there  exists  a vector 

/ II 

k such  that  c = k X X . 

(iii)  Y is  estimable  if  and  only  if  c = c H,  where 
H = (x'x)"  x'x. 

It  should  be  noted  that  in  practice  it  is  not  a trivial 
matter  to  check  for  estlmability  due  to  complicated  nature  of 
the  design  matrix  X.  This  is  why  the  experimenter  is  well 
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i 

advised  to  specify  his  set  of  linear  parametric  functions  of  i 

interest  and  try  to  collect  his  data  (i.e.,  chooses  his  de- 
sign matrix  X)  which  guarantees  the  estimability  of  his 
functions  of  interest. 

Definition  We  say  the  design  X is  connected  for  Y if 
'y  is  estimable  under  X.  Otherwise,  X is  said  to  be  dis- 
connected for  Y. 

I 

Estimation,  Test  and  Confidence  Interval  for  Y.  j 

If  Y is  estimable  under  X then  it's  well  known  that 

the  bes-t  linear  unbiased  estimator  of  Y is  given  by 

^ /A  / / _ / 

Y = c 3 = c (X  X)  X Y, 

with 

var  '1^  = = c (X  X)”  ca^ . 

V 

An  unbiased  estimator  of  is  c (X  X)  cc  where 

v I 

y'(I-X(x'x)“X^)Y.  j 

Therefore  * I 

^ ~ N(»,  c\).  I 

T I 

We  know  that  j 

Ql  = ^ Y'(l-X(x'xr  X')Y  = §2  ^ ^2  | 


T is  a linear  form  in  Y and  thus  Y and  Q,  are  indepen- 


IJIH  l,n.  nwpi^ipwpi 


1 
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dent  If 


c'(x'x)-  x'[i-x(x'x)~  x']  = c'(x'x)"  x'  - c'(x'x)“  x'x  (x'x)"  x'  = 0.'^ 

I I 

Now  since  Y is  estimable  thus  c can  be  expressed  as  i 

• II  ! 

t X.  Therefore,  by  substituting  t X for  c in  the  above  ‘ 

I 

expression  we  obtain: 


/ / 


' / ' X - ' / ^ ' 


= t X(X  X)“  X - t X(X  X)"  X X(x  X)"  X 


/ / 


= t X(X  X)"  X - t X(X  X)"  X = 0, 


Since  T r-  N(Y,  cta),  this  imples  that 


A 

Y - Y 


Y - Y 


Thus 


^(c'(x'x)-  c)l/2 


(^-Y)/g(c'(X^Xr 

[Q^/(n-r)]l/^ 

Y-  Y ^t(n-r). 


N(0,  1)  . 


J [c'(x'x)- 


A 

or  equivalently 


^ Y - Y ^ ^ F(l,  n-1) 


A 

This  statistic  can  be  used  for  testing  hypothesis  of  the  form 
Ho:  Y = m.  This  statistic  can  also  be  used  for  constructing 


■J 
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confidence  intervals  for  Y. 


or 

P ^ Y ^ + aA(F^(l,n-l))^/^]  = 1 - a. 

Suppose  we  have  a set  of  linear  paramteric  functions 
Yg > • • • wish  to  construct  1 - a simultaneous 

confidence  intervals  for  these  t linear  parametric  functions ^ 
The  above  confidence  intervals  give  1 -a  confidence  inter- 
vals for  individual  Y's.  Scheffe's  S-raethod  answers  this 
problem. 

But  first  we  need  the  following  definitions. 

Definition  3.  Two  linear  parametric  functions  Y^  and  Y^ 
are  said  to  be  algebraically  independent  if  their  correspond- 
ing coefficient  vectors  c^  and  c^  are  independent. 

Definition  4.  By  a q-dimensional  subspace  L of  linear  es- 
timable functions  under  the  design  matrix  X we  mean  the 

subspace  generated  by  the  coefficient  vectors  of  q indepen- 

/ 

dent  lieanr  estimable  functions  under  X.  We  say  Y = c p e L 
if  c^  e L. 
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The  S-method  of  multiple  comparison  is  based  on 

Theorem  1..  Under  Q the  probability  is  1-a  that  all 
estimable  functions  Y in  a given  q-dimensional  space  L 
simultaneously  satisfy 


(1)  ^ ^ '1'  ^ ^ , 

where  S = [q  F (q,n-r)]^^^. 

One  can  rewrite  (l)  in  the  following  form 


Ft  l^-^l  < [q  F^t^q,  n-l)]^/^  for  all  Y e L]  = 1 - a, 
or  2 

p^(^-Y)  /c'jX  X)~  c ^ F (q,n-l)  for  all  Y e Lj  = 1 - a. 

A2 

a 

A 'A  / 

Since  Y=cp,  Y=cp  therefore  it  suffices  to  prove  that 

the  maximum  value  of  (c  ^ - c 3)/c  (X  X)”  c for  all  nonzero 

2 2 

c € L is  distributed  as  a X (q);  and  this  maximum  is  inde- 

^2  2 2 
pendent  of  (n-l)a  which  is  distributed  as  a X (n-r). 

To  prove  this  we  need  the  following  lemmas. 


Lemma  1.  Let  A be  a symmetric  matrix  of  order  n.  The 

/ / 

maximum  value  of  z Az/z  z over  all  nonzero  z e is 
the  largest  eigenvalue  of  A,  and  this  maximum  is  attained 
when  z is  any  eigenvector  of  A corresponding  to  the  root 

t. 


I 


|3 


ii 


I 


LL 
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Froof . First  we  show  that  the  following  two  problems  are 
equivalent . 

(i) 


maximize  5'  ^ (Fi)  maximize  z Az  . 

z40  z'z  z' z=l 


Let 


max 

/ . / 

Z -fiz  = tn  and  max  7.  A 7 = 

m2 

z40 

z'z  ^ z'z=l 

and 

suppose 

m^  is  attained  for  z = z^ 

and  mg  is  attained 

for 

Z = Zg, 

i .e . , 

/ 

^ I 

1 1 = m^  and  z^Azo  = 

Z-Zi 


igAzg  - m^ 


Let 


"l  = 


1 _ ' - 

Zt , then  z^  z,  = 1, 

/_  1 11 


thus 


_ / _ 

z^Az^ 


ZiZi 


z- Az^ 

= m^, 


Therefore,  m^  < m^ . Also  since  ZgZ^  = 1,  Zg  4 0 and 


1 - "'s 


^2^2 


thus  m^  ^ m^.  Hence  m^  = m2 > 


Note ; Since  max  z Az/z  z is  equivalent  to  max  z Az  and 
z40  z^z=l 


/ / 

z Az  is  a continuous  function  of  z and  lz:zeE,zz=lJ 


n' 


is  a compact  set,  it  follows  that  max  z Az/z  z exists 

z40 


•a 


i 


J 


We  shall  now  present  two  methods  of  finishing  the  proof  of 
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Lemma  1 . 


Method  1 (Lagrange  multiplier  method). 

/ / 

We  want  to  maximize  z Az  subject  to  the  constraint  z z = 1. 


f(z,X)  = z Az  - X.(z  z-l) 

S'f* 

where  \ is  a Lagrange  multiplier.  Since  — "2Az-2  Xz=0. 

oZ 

this  imples  that  Az  = Xz  so  X must  be  an  eigenvalue  of  A. 

On  the  other  hand,  if  X is  an  eigenvalue  of  A,  then 


so  that 


' / r f 

z Az  = z (Az)  = z (Xz)  = Xz  z. 


/ / 

max  z Az  = max  Xz  z = max  X 
/ / 

Z Z=1  Z Z=1 


where  X is  an  eigenvalue  of  A.  Therefore 

/ / 

max  z Az  = max  z Az  = max  \ - Z . 
z4=0  z^z=l 


Now  suppose  V is  any  eigenvector  corresponding  to  Z,  then 
V Av  - V ( Av)  =v  /v_  Xvv_ 

I I • I 

vv  vv  vv  vv 

so  that  the  maximum  value  is  attained  for  any  eigenvector 
associated  with  Z. 

Method  2.  Since  A is  a real  symmetric  matrix  of  order  n, 

/ 

there  exists  an  orthogonal  matrix  P such  that  P AP  = A 
where  A is  a diagonal  matrix  with  the  (real)  eigenvalues 
of  A on  the  diagonal.  Let  P.  denote  the  j-th  column 
of  P,  i.e.,  P = (P, :P„ ; . . . ;P  ) , then  P.  is  an  eigenvalue 

* La  • • n yj 
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of  A corresponding  to  X satisfying 

J 


p^p  ^ fO  when  j 4 
J k ll  when  J = k 


It  follows  that 


A = pap'  = + ^pPpPp  +...+  ^ P p1 

111  L c c nnn 


I = pp'  . p^p'  + p^p^  p^p'  . 

The  set  {P^,P„,...,P  } is  an  orthonormal  V -sis  for  E . 

J.  ^ n 

Let  z € E^,  then  z = Pw,  where  w = (w^,'- ^ • ^w^)  and 
w^'s  are  the  coordinates  of  the  vector  z witn  respect  to 
the  basis  [P^^Pg^ • • • ^P^l • Therefore 


z'az  _ (Pw) 'pap' (Pw)  _ w'p'pAp'pw 


z ^ z ( Pw  ) ' ( Pw  ) 


w P Pw 


' X,W^  + X„w^  + . . . + X w^ 

_ w Aw  _ 11  2 2 n n 


w w + Wg  + . . .+  w^ 

/ / 

So  maximizing  z Az/z  z for  z 4 Oj  2;  e E^  is  equivalent 


to  maximizing 


XjW^  + XgWg  f . . .+  X^w^  w € E^ 

w|  + w|  + . . .+  w^ 


and  w 4 0 (since  z 4 0) • 

Suppose  max  ( X^, Xg > • • • ^ X^)  = then 


^l”l  ^2'"2  ■*■•••+  ^n”n  - ^ i=l 

2 2 2 
X,w,  + X,w„  + ...+  X i 

1 1 c c n n ^ _ 

2 2-2  " 

Wi  + Wg  + . . .+ 


^ 2 
A.S.wf 
1=1  1 

^ 2 
F.  w. 

i=l  ^ 


If  jI  = X . then  let  w = 0 if  i =}=  J w . 4 0*  Then 

J 1 .1 


2 2 2 
X-Wt  + X„w„  +...+  X 
11  c c n n 


Wt  + w?  + . . .+  vj2 

12  n 


so  i is  attainable  and  the  maximum  is  attained  for 
w = {O, . . . ^0,w  ,0, . . . ,0}  or  equivalently  for 

cl 

z = Pw  = P.W.,  w.  4 Oj  i.e.,  for  any  eigenvalue  of  A 
JO  J 

corresponding  to  I the  largest  eigenvalue  of  A. 


Lemma  £.  Let  A be  a symmetric  matrix  of  order  n and 

let  B be  any  positive  symmetric  matrix  of  order  n.  The 

/ / 

maximum  value  of  z Az/z  Bz  over  all  nonzero  z e E is  -t, 
the  largest  eigenvalue  of  B“^A,  and  this  maximum  is  attained 
for  any  eigenvector  of  B"*^A  corresponding  to  the  root  I . 

Proof . First  let  us  solve  the  following  problem.  Maximize 
z Az  subject  to  z Bz  = 1.  Using  the  method  of  Lagrange's 
multiplier  let 

f(z,X)  = z Az  - X(z  Bz-l) 

where,  X is  a Lagrange  multiplier.  A necessary  condition  is 
3f/az  = 2Az  - 2XBZ  = C 


7 
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r 


Az  = XBz  = B\z  =»  B”  Az  = \z . 

So  that  X is  an  eigenvalue  of  B~^  and  z is  the  corres- 
ponding eigenvector.  Therefore, 


max 

/ 

z Bz=l 

/ / 

z Az  = max  z (Az)  = 

z Bz=l 

max  z (XBz) 

/ 

z Bz=l 

— 

max 

/ 

Xz  Bz  = max  X = i . 

z Bz=l 

Now  we 

shall 

show  that  the  following 

two  problems  are  equiva 

lent : 

(i) 

max 

zeE 

n 

z^O 

/ 

z Az  . ( ii)  rpax 

z Bz  z Bz=l 

/ 

z Az 

Proof  of  this  is  very  much  like  the  counterpart  proof  given 
in  the  beginning  of  Lemma  1. 

Let 

max  = m,  and  max  z Az  = m2 

z+0  z^Bz  z'bz=1 

/ » 

and  suppose  is  attained  for  z = v^  and  1^2  is  attained 

for  z = V2 , i.e.. 


m^  and  ^2^^ 2 “ 


Let 


^1  = 


(v^Bv^) 


1/2 


-11- 


r 

i 


r 

I 


\ 

r 

t 

1 

i- 


then 


also 


_/  - v-Bv, 
''1®''-!=  = 1, 
ViBVi 


v[av^  _ v'av^  _ ^ 

— ■ —7 ^ • 


v^BVi 


Therefore,  in,  < nu, . On  the  other  hand,  since  v^Bv_  = 1 

1-^^  Ik  d 2 


and 


since  B is  positive  definite 


4 


''2®''2 


v'aVj 


= 


thus 


Hence  . 


We  shall  shortly  give  a generalization  of  the  preceeding 
results.  Let  L be  a vector  subspace  of  of  dimension 

q.  Let  the  columns  of  C = [ c^,C2 , . . . ,0^]  be  a basis  for  L. 

/ / 

Lemma  3.  The  maximum  value  of  z Az/z  z over  all  nonzero 
z e L is  I,  the  largest  eigenvalue  of  CC'*’A,  and  is  attained 
when  z is  any  eigenvector  of  CC'*’A  corresponding  to  the 
root  t,  where  C'*’  is  the  Moore-Penrose  generalized  inverse 
of  C. 

Proof . C is  an  n x q matrix  of  rank  q thus  (C  C) 
exists  and  one  can  check  that 

= (cc)  c , 

is  the  Moore-Penrose  generalized  inverse  of  C.  Also  note 
that  matrices  CD  and  DC  have  the  same  eigenvalues.  Now 


let  V e Eq,  then  z = Cv  e L,  hence 


max  14s  > max  iP-'') 

zeLj  z40  z z v40  (Cv)  (Cv) 


/ / 

max  V (C  AC)v 


v40  ^'(c'c)v  ■ 


Since  C G is  positive  definite  we  can  use  Lemma  2 and  con- 
c lude  that 

/ / 

max  V (C  AC )v 
veEq,v40  v^(C  C)v 

= largest  eigenvalue  of  (C  C)~^  C AC 

= largest  eigenvalue  of  C'*’AC 

= largest  eigenvalue  of  C(C‘*’A)  (see  the  re- 
mark about 
CD  and  DC ) 

= 1. 

To  obtain  the  inequality  in  the  other  direction,  let  z e L. 
This  implies  that  tnere  exists  a v e E^  such  that  z = Cv. 
Thus 


max  z Az  ^ 


zeL,  z40 


max  (Cv)  A{Cv] 


''+°  (Cv)'(Ov) 


= largest  eigenvalue  of  CC'*'A  = t. 


Now  let  z be  any  eigenvalue  of  CC'*’A  corresponding  to  the 

root  -t . Then  CC  Az  = Az  which  imples  z CC  Az  = Az  z 

/ / 

which  in  turn  implies  z Az/z  z = A because  z e L implies 


that  z CC  = z . 
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Corollary  1.  Let  H be  any  q x n matrix  of  rank  q.  Then 

/ / 

the  maximum  value  of  z Az/z  z over  all  nonzero  z in  E 

n 

satisfying  Hz  = 0 is  i- , the  largest  eigenvalue  of 
and  is  obtained  when  z is  any  eigenvector  of  (l-H'^H)A  cor- 
responding to  the  root  1.  Here  H'*'  denotes  the  Moore-Penrose 
generalized  inverse  of  H. 

Proof . Hz  = 0 if  and  only  if  z belongs  to  the  column 
space  of  I-H'*’H.  This  is  seen  as  follows: 

If  Hz  = 0,  then  z is  in  the  coliunn  space  of  I-H'*'H,  i.e., 
there  exists  a w such  that  (I-H‘*’H)w  = z.  Set  w = z then 
(I-H'*'H)z  = z - H'^Hz  = z - H^(Hz)  = z - H^(0)  = z.  On  the  other 
hand,  if  z is  in  the  column  space  of  I-H'*’H  then 
Hz  = H(l-H'*'H)w  = Hw  - HH‘*'Hw  = Hw  - Hw  = 0 . Now  the  proof 
follows  from  Lemma  3. 

Lemma  4.  Let  B be  a positive  definite  matrix  of  order  n. 

/ / 

Then  the  maximum  value  of  z Az/z  Bz  over  all  nonzero  z in 

L is  I,  the  largest  eigenvalue  of  C(C  BC ) C A,  and  is 

/ / + » 

attained  when  z is  any  eigenvector  of  C(C  BC)  C A corres- 
ponding to  the  root  h . 

Proof.  An  argument  similar  to  the  one  used  in  the  proof  of 
Lemma  3 gives  us 

max 

zeL,  z^O 


I 

I 

I 

\ 

) 

,1 

ii 


I 


By  an  earlier  result  one  gets 


• • • JL  • 

max  V C ACv  _ largest  eivenvalue  of  [(C  BC)'*’  C AC] 

''+°  v'c'bCv 

i 1 / 

= largest  eigenvalue  of  [C(C  BC)*  C A]  = i 
(recall  the  argument  about  CD  and  DC). 

Now  let  z be  any  eigenvector  of  C(C  BC)'*'  C A corresponding 

to  the  root  1.  Then 

/ 1 / 

C(C  BC)  C Az  = -tz  which  implies  that 

z BC(C  BC)"'"  C Az  = Iz  = Iz  Bz,  which  implies  that 

/ / 

z Az/z  Bz  = /,  since  z e L implies  that 
z BC(C  BC)^  C = z . 

This  latter  claim  is  seen  as  follows.  Since  z e L then 
there  exists  a w such  that  Cw=z,  i.e.,  z is  a linear 
combination  of  the  columns  of  C which  generate  L.  Then 

z BC(C  BC)^  C = w C BC(C  BC)^  C = w C = (Cw)  = z . 

Lemma  5.  The  Moore-Penrose  of  X is  given  by  (X  X)"^X  where 
A"^  denotes  the  Moore-Penrose  of  the  matrix  A. 

Proof.  By  definition  K is  the  Moore-Penrose  generalized 

inverse  of  A if  AKA  = A,  KAK  = K,  (AK) ' = AK  and 

/ 

(KA)  = KA.  Therefore,  we  shall  check  these  four  conditions 
for  X'*’ . In  what  follows  we  use  the  following  well  known 
facts : 
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F^:  X(X  X)”X  is  symmetric  and  X(X  X)^X  = X(X  X)"X  where 
(X  X)”  is  any  generalized  inverse  of  (X  X). 

Fg:  X(x'x)"x'x  = X and  x'x(x'x)“x'  = x'. 

(i)  xx"^x  = x(x'x)"^x'x  = x(x'x)"x'x  = X, 

(ii)  x'^'xx'’'  = (x*x)'^x*x(x'x)‘^x"  = (x'x)'*'x'x(x'x)"x"  = (x'x)'^x'. 

(iii)  (XX"^)'  = (X(x'x)"^x') ' = (X(x'x)"x')'  = X(x'x)"  = X(x'x)'^x' 


= XX'*', 


\ + 


(iv)  (X“^X)  = ((X  X)  X X)  = (X  X)  since  (X  X)  is  the 

Moore-Penrose  inverse  of  XX  and  thus  (XX)  XX 
is  symmetric. 

Lemma  6.  The  Moore-Penrose  generalized  inverse  of  X is  (X  ) 
Proof . 

(i)  x'(X'*')'x'  = [XX'*'X]'  = [X]'  = x', 

(ii)  (X'*')'x'(X'^) ' = [X^’^XX'*']'  = [X'*']', 

(iii)  [x^x-*-)']'  = [(x'*'x)']'  = [x'*'x]'  = x'(x-^)', 

(iv)  [(X'*')'x']'  = [(XX'*')']'  = [XX"^]'  = (x'*')'x'. 


Lemma  2*  Moore-Penrose  of  X,  then  XX"*'(X'*') 

= (X*")'. 

Proof.  From  Lemma  5 X'*^  = (X  X)'*^X  . Thus 

xx'*'(x'*')'  = x(x'x)'*'x'  [ (x'x)'‘'x'] ' 

= x(x'x)'*'  [x'x)‘''x'x] ' 

= x(x'x)'*'  [x'x[  (x'x)'*'] '] 

= X(x'x)'*'  [x'x[x'x) ']'*■]  by  Lemma  6 
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= X(X  X)^  [X  X(X  X"*"]  by  a property  of  Moore- 

Penrose 


= X(X  X)'*’  generalized  inverse 

= [[(X  X)"^]  X ] = [[(X  X)  rx  ] 

= [(x'x)‘^x']'  = (X'^)'. 

Lemma  8.  The  set  of  nonzero  eigenvalues  of  DC  coincides 
with  the  set  of  nonzero  eigenvalues  of  CD. 

Proof.  Let  and  A^  be  the  set  of  nonzero  eigenvalues 

of  DC  and  CD  respectively. 

If  (DC)x  = Xx  =»  C(DC)x  = XCx 
=>  CD(Cx)  = X(Cx)  = CDy  = Xy, 

so  if  X is  an  eigen-'-^alue  of  DC  it  is  an  eigenvalue 
of  CD, 

=»  A^cAg . Similarly,  Ag'^A^.  Thus  ~ ^2* 


Proof  of  Theorem  1. 


( 1) . max  (c 
ceL  7^ 


a'xl-a'x6)^ 


o'(x'xro"  a'x(x'x)-x'a 


Reason,  c P is  estimable  =*  a an  a such  that 

/ / / III 

c = a X.  But  c e L =*  c = D t^c^,  c^  = L^X,  c^’s 


were  chosen 


c = E t^b^X  = [Et^bj^JX  = a X so  a 


is  a linear  combination  of  b^*s.  But  from 

II  II  / 

Cf  = b^X  =»  [Et^b^JX  = a X so  a is  a combination 
of  b^'s.  But  from  c^  = b^X  ^ = c^^^  or  = (X 


(a  XX  Y-a  Xp)  . Reason.  See  Lemmas  5 and  6. 

a'xx"^a 

(3)  (a'xx"^Y-a'xf3)^  ..  (a'y-a'xs)^  a ' ( Y-Xp  ) ( Y-Xp ) 'a 

> ~ / “ / • 

a XX  a a a a a 

Reason.  Since  a e ■£[  (X  )"^C]  ^ a e column  space  of 
(X  i.e.,  3 an  f such  that  a = (X  )’*’f  = (X"^)  f 

=»  a = f X . Thus  a XX^  = f X XX  = f X^  = a . 

(4)  From  (l)  and  (3) 

max  ( c P-c  S _ max  a (Y-Xp)^Y-Xp)  a 

c (X  X)"c  ae.C[  (x')'^C]  a a 

= max  & k&  , A=  (Y-Xp) (Y-Xp) ' 

aer[  (x')'^C]  a a 

= largest  eigenvalue  of  [ (X  ' )'^C][  (X  ' )'*'C]‘^(  Y-Xp  ) (Y-Xp ) ' 
by  Lemma  3 . 

= largest  eigenvalue  of  (Y-Xp)  [(X  )'*^C][(X  )'^C  ]"^(  Y-Xp ) 
by  Lemma  8,  but  this  is  a scalar, 

= (Y-Xp) '[  (x')'^C][  (x')'‘’C]'^(Y-Xp)  = which  is  a 
quadratic  in  (Y-Xp)  ~ N(0,  o^I) . 

The  claim  is  that  Q ~ This  will  be  the  case 

if  we  prove  that  [(X  )’*’C][(X  )^C)"'’  is  idempotent  and 
its  rank  is  q.  The  Idempotency  is  obvious  since  in 


_ ^ 
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general  = BB’^BB'*'  = BB"*^  . We  shall  now  show 

that  rank  [(X  )‘^C][(X  = q.  This  can  be  seen 

as  follows : 

r[(x')+C][(x')+C]'^  ^ r[(x')‘^C]  ^ r[C]  = q, 
on  the  other  hand, 

r[{x'x)‘^C][(x')+C]‘^  > r[(x')‘^C][(x')+C]'^[(x')'^C] 

^ r[(x')'*'C]  = r[(X‘^)'c]  = r[  [ (x'x)'^X ' ] 'c]  = r[X[  (x'x)"^] 'c] 

> r[x'x[  (x'x)"^] 'c]  = r[x'x[  (x'x)"^] 'x'k]  since  C = x'k 

= r[[X(x'x)‘^x'x]'K]  = r[[X)x'x)"x'x)'K]  = r[x'x(x'x)"x'K] 

= r[x'K]  = r[G]  = q. 

The  proof  of  Theorem  1 will  be  complete  if  we  show  that 
and  Qg  are  independent  where, 

Qg  = (n-r)  = y' ( I-X)X 'x) “X ' )Y 

= (Y-Xp) '(l-X(x'x)"x')(Y-XB) . 

It  is  sufficient  to  prove  that 

[i-x(x'x)'x'][(x')'^c][(x')'^c]^  = [i-x(x'x)"x'][(x'^)'c][(x+)'c]'^  = 

by  Lemma  6 

LHS  = [ (X”^) 'c][  (X'*') 'c]"*"  - X(x'x)"x'(X'*‘) 'c[  (X"^) 'C]"^ 

= [ (X'*') 'c][  (x"^) 'c]*^  - x(x'x)'^x'(x'^) 'c[  (C"^) 'c]'* 

= [ (X"^) 'c][  (X'*') 'c]"^  - XX‘^(X'^) 'c[  (X'^) 'c]"^ 

= [(X+)'c][(X'^)'c]  - (X‘*')'c[(X+)'c]  = 0. 


The  relation  of  the  S-method  or  S-intervals  for  '1'  in  L 
and  the  standard  F-test  of  the  hypothesis 

«o=  h = ’z  = = 0 

is  stated  in 

Theorem  2.  Under  fl  the  a- level  F-test  of  Hq  will 
accept  Hq  if  and  only  if  for  all  Y in  L the  intervals 
(l)  in  Theorem  1 cover  zero. 
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